Detection of 33.8 ms X-ray pulsations in SAX J0635+0533 
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ABSTRACT 

We revisited the BeppoSAX observation of SAX J0635+0533, suggested as a 
possible counterpart to the gamma-ray source 2EG J0635+0521. We have discovered 
a 33.8 ms pulsation from the source and derived an improved position, consistent 
with the location of the Be star proposed as a binary companion. We interpret 
the periodicity as the spin period of a neutron star in a binary system with a Be 
companion. 

Subject headings: X-rays: stars — pulsars: general — gamma rays: observations — 
stars: individual (2EG J0635+0521, SAX J0635+0533) — stars: neutron 
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1. Introduction 

The X-ray source SAX 10635+0533 was discovered by Kaaret et al. (1999) thanks to a 
BeppoSAX observation within the error box of the unidentified Galactic gamma-ray source 



2EG J0635-I-0521 ( ^Thomson et al. 1995| ), a candidate gamma-ray pulsar as suggested by its hard 



gamma-ray spectrum ( [Merck et al. 1996| ). The X-ray source is characterized by quite hard X-ray 



emission detected up to 40 keV ( pCaaret et al. 1999] ). Its energy spectrum is consistent with a 
power-law model with a photon index of 1.5, an absorption column density of 2.0 x 10^^ cm"^, 
and a flux of 1.2 x 10"^' ergcm"^s"^ in the 2-10 keV energy band. A search for pulsed emission 
over a period range from 0.030 s to 1000 s did not detect any pulsed signal. Due to the large error 
box of the gamma-ray source, the identification of SAX J0635-I-0533 with 2EG J0635-I-0521 is not 
definitive: such an identification could only be made through pulsed detection in both X-ray and 
gamma-ray emission or a much improved gamma-ray position. 



Follow up optical observations ( [Kaaret et al. 1999| ) suggest as a counterpart of 



SAX J0635-I-0533 a Be star with a V-magnitude of 12.8, located within the 1' X-ray 
source error box. The estimated distance is in the range 2.5-5 kpc. The total Galactic 21 cm 
column density along the optical counterpart direction is 7 x 10^^ cm"^ ( [Stark et al. 1992[ ), in 



agreement with an estimation derived from the extinction of the optical spectrum. Kaaret et al. 
(1999) assert that the larger column density for SAX J0635-I-0533, estimated from the X-ray 
spectrum, implies the presence of circumstellar gas around the X-ray source. Moreover, taking 
into account the positional consistency between the X-ray and the gamma-ray sources as well as 
the hard X-spectrum, the strong X-ray absorption, and the optical association with the Be star, 
they suggest that SAX J0635-I-0533 is an X-ray binary emitting gamma-rays. 

We revisited the BeppoSAX observation of SAX J0635-I-0533, available from the public 
archive (obs.code #30326001). In this letter we present new imaging and timing results. Our 
analysis has revealed a 33.8 ms pulsation of the X-ray source. §2 describes the observation and 
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data reduction. In §3, we report the data analysis procedures and results. We conclude with a 
discussion of the results in §4. 



2. Observation and Data Reduction 

The field including the serendipitous source SAX J0635+0533 was observed on 23-24 
October 1997 by the Narrow Field Instruments (NFIs) on board BeppoSAX dBoella et al. 1997a[ ). 



In our analysis, we use only data coming from the two NFI imaging instruments, namely the Low 
Energy Concentrator Spectrometer (LEGS) operating in the energy range 0.1-10 keV (P armar et 



al. 1997) and the Medium Energy Concentrator Spectrometer (MECS) operatmg m the energy 



range 1.3-10 keV ([Boella et al. 1997B[). At the observation epoch only two of the three MECS 



detector units were operating. 

LECS and MECS raw data have been reduced in cleaned photon list files by using the 
SAXDAS v.2.0 software package and adopting standard selection criteria^. The source was 
located ~ 3' off-axis. The total duration of the observation is 72714 s and the net exposure is 
14043 s and 34597 s for LECS and MECS, respectively. The different exposures are due to the 
fact that the LECS was only operated during satellite nighttime. 



3. Data Analysis 

In our analysis, we include only events within a 3.7' radius around the source position; this 
radius optimizes the signal-to-noise ratio for SAX J0635-I-0533 and contains ~80% of the source 
signal for both LECS and MECS. Due to the fact that the source is strongly absorbed, we further 
improve the signal-to-noise ratio by using the energy range 1.8-10 keV for both the instruments. 



' http ://ww w. sdc. asi.it/software/saxdas 
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The resulting number of selected events is then 479 for LEGS and 3 658 for MECS, respectively. 
We estimate that less than 1% of these events is due to the diffuse and instrumental background. 

3.1. Spatial Analysis 

The position of SAX J0635+0533 has been re-evaluated by considering MECS data 
accumulated under specific constraints of the spacecraft; in particular, only data referring to the Z 
star tracker (the one aligned with the NFIs) in use have been taken into account. This allows us to 
reduce the attitude reconstruction error to ~ 30" (90% confidence level)0. 

After suitable smoothing of the screened data and then by using a Gaussian centering 
procedure, the position has been found at i?A = 06''35"'18', DEC = 05°33'11" (J2000) with no 
significant statistical error. The error box is only determined by the attitude reconstruction error. 
The Be star (Kaaret et al. 1999) is within the 30" of the SAX J0635+0533 error box and 6.8" 
distant from its center. 

EDITOR: PLACE HGURE |l] HERE. 

Fig. 1 shows a finding chart for SAX 10635+0533 from the Digitized Sky Survey in a field 
of 4' X 4'. Within the reduced X-ray source error circle, only two of the seven stars referred to in 
Kaaret et al. (1999) are now present. The Be star is the brightest one which is positioned close to 
the center of the X-ray source error circle. 



^http : //www. sdc . asi. it/software/cookbook/attitude .html 
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3.2. Temporal Analysis 

The SAX J0635+0533 light-curve (1000 s bin size) for the MECS is shown in Fig. 2a 
(gaps are present due to non-observing time intervals during South Atlantic Anomaly and Earth 
occultation). Fig. 2a shows that the emission of SAX J0635+0533 is variable up to a factor of 
10. Fig. 2b shows the hardness ratio between 4-10 keV and 1.8^ keV bands. The relative 
contributions of these two energy bands are not constant, but no strong correlation with respect to 
the total intensity is evident. 

EDITOR: PLACE FIGURE |HERE. 

In order to search for periodicity, the arrival times of all selected events have been converted 
to th e Solar System Baricentric Frame, using the BARICONV codefj. The Zf test (B uccheri et al. 
1983) on the fundamental harmonics with the maximum resolution {6f = 1/72714 Hz) applied to 
the MECS baricentered arrival times does not reveal significant deviations from a statistically flat 
distribution up to 50 Hz. 

If SAX J0635+0533 is a binary pulsar of rotational spin period and orbital period Pg, the 
observed Ps is modulated by the orbital motion. Thus, a direct search for a coherent oscillation at 
Ps can be successful only if the modulation amplitude is small over the time interval AT in which 
the search is performed. This condition is satisfied if AT <C /'o- To reduce the effect of a possible 
orbital motion in the periodicity search, we divide the whole data span into M subintervals, 
calculating the Zj statistics for each trial period in each subinterval, and then adding together the 
M statistics for each trial period. This procedure results in a less noisy spectrum. The sum of the 
M separate Z^ statistics results in a statistics with 2M degrees of freedom ( |Bendat & Piersol 1971| ) 
which we refer to as the Zf (z/ = 2M) statistics. Because the power depends on the square of the 

^ |http://www.sdc. asi.it/software/saxdas/baryconv.html 
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pulsed signal, its strength decreases with M and only sufficiently strong signals can be detected. 
We selected time slices corresponding to intervals of continuous observation taken between two 
Earth occultation periods. The total number of these slices is M = 13, each one lasting ~ 3300 s. 
We adopted a frequency step 5f = 2x 10"^ Hz, spanning 50 Hz of search range with 250000 trial 
frequencies. 

EDITOR: PLACE FIGURE ^HERE. 

Fig. 3 shows the power spectrum obtained with the MFCS and LEGS data in the Z^iu = 26) 
statistics as a function of frequency, where an evident excess appears at /o = 29.5364 ±0.0001 
Hz. The value of Z^(z/ = 26) at this frequency is equal to 99.6. Because the Z^(z/ = 26) follows the 
statistics with 26 degrees of freedom, the single trial chance occurrence probability to have an 
excess greater than 99 is 2 x 10"^", as shown in Fig. 4. Taking into account the number of trial 
frequencies used, the probability is 5 x 10"^, corresponding to 4 standard deviations of the Gauss 
statistics. When we use only MEGS data, the Zj{u = 26) value decreases to 91 in agreement with 
the reduction of the source counts. We tested also the power at half and double values of the 
detected frequency /q: no signal is present at / = /o/2 or at / = 2/o. In order to check for the 
persistence of the periodicity in the whole observation, we also binned the maximum resolution 
spectrum of the entire data set by a factor 13. We obtain a high {u = 26) > 104, at frequency = 
29.5364 Hz, indicating that the previous high power value does not come from one or few selected 
time intervals. 

EDITOR: PLAGE FIGURE |HERE. 

From the Zj (z/ = 26) value we can estimate the pulsed fraction. For Np pulsed counts over 
Nt total counts, Z^(z/ = 26) = 2aN^/Nt + v, where a is a shape constant. In our case, a = 0.25 
(sinusoidal shape), and the pulsed fraction is then about 0.2. 
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We also performed a periodicity search over the full observation including the first and 
second derivative of the frequency, / and / respectively, at the maximum resolution step. The 
search was performed within the interval A/ = 2 x 10"^ Hz centered on the detected value /q. 
The value range for / and / were chosen consequently. We obtained a maximum in the power 
spectrum (Z2(zy = 2) = 52) for / = 29.53643 ± 0.00001 Hz, / = (-3. 1 ± 0.2) x lO'^ Hz s"^ and 
/ = (1.1±0.1) X 10"'"^ Hz s"^, where the errors refer to the parameter resolution. Assuming a 
circular orbit, the second order polynomial expansion of the frequency vs. time relation is locally 
compatible with a orbital period ~ 18 days and semi-major projected axis apSinz ~ 63 s It. 
However, we stress that the interval of observation is much shorter than the derived orbital period, 
and that the parabolic fit may not be a good representation if the orbit is eccentric. Thus, these 
results should be taken only as a possible indication of the orbital parameters and not as a firm 
detection of orbital motion. 

Figure 5 shows the pulse profile folding the data taking into account the frequency derivative 
terms. Note that if the orbit is eccentric, the real pulse profile could be narrower. 

EDITOR: PLACE FIGURE | HERE. 



4. Discussion 

The coherence of the detected periodicity is high, Q = f / 5f ^ 10^. This value is much 
greater than those observed in quasiperiodic oscillations (QPOs) often detected in the X-ray 
emission of X-ray binaries ( [van der Klis et al. 1996| ). The high coherence induces us to interpret 



this periodicity as a neutron star spin period. 

The association between SAX J0635-I-0533 and the Be star has been reinforced thanks to a 
reduced error circle of the X-ray position coordinates. Timing analysis results indicate that the 
neutron star could orbit around a companion star. The tentative orbital parameters are consistent, 
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for orbit inclination less than 25°, with a primary mass greater than 10 Mq, as expected for a 
Be-star. 

The X-ray emission may be powered either by accretion or by spin-down of the neutron star. 
We consider these possibilities in turn. 

The SAX J0635-I-0533 system may consist of a rotation-powered pulsar orbiting the Be star. 
In this case, the X-ray emission could be magneto spheric emission similar to the power-law 
component in the X-ray emission of known X-ray /gamma-ray pulsars. The pulsation frequency 
we have found and the X-ray luminosity of the source are similar to those of the known 
X-ray/gamma-ray pulsars. The high X-ray variability of SAX J0635-I-0533 on time scales of 
1000 s is unlike the steady X-ray emission seen from isolated pulsars. However, it could be 
produced by variable X-ray absorption caused by matter in the binary system or a wind from the 
Be star. 

An alternative, but still rotation-powered, scenario is that SAX J0635-I-0533 is similar to the 
Be radio pulsar PSR J 1259-63. These two sources have similar spin frequencies and similar X-ray 
spectra ( [Nicastro et al. 199^ ). In this case, the X-ray emission of SAX J0635-I-0533 should arise 
from a shock interaction of the energetic particles from the pulsar with the wind from the Be star. 
However, the upper bound, 8%, on the X-ray pulsed fraction from PSR J 125 9-63 in the 2-10 keV 
band (Kaspi et al. 1995) is well below the value estimate for SAX J0635-I-0533 in the same energy 
band. 

The X-ray emission from SAX J0635-I-0533 may be powered by accretion. Strong X-ray 
variability would naturally occur in such a system. Following this interpretation we can infer the 
magnetic field strength of the neutron star. For accretion to proceed, the centrifugal force on the 
accreting matter co-rotating in the magnetosphere must be less than the local gravitational force 
dlllarionov & Sunyaev 1975] ; [Stella, White, & Rosner 1986| ). Assuming a bolometric luminosity 



of 1.2 X 10 ergs (0.1-40 keV) estimated from spectra results given in Kaaret et al. (1999) for 
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a 5 kpc distance, and a neutron star mass of 1.4Mq and radius of 10 km, we can set an upper 
limit on the magnetic field strength of 2 x 10^ G. This is a factor 10^ lower than those measured in 
typical accreting X-ray pulsars, but similar to the fields inferred for the 2.49 ms low-mass X-ray 
binary SAX J1808.4-3658 dWijnnands & van der Klis 1998D and for millisecond radio pulsars. 



The X-ray luminosity of SAX J0635-I-0533 is a factor of 10 below that of most Be/X-ray binaries 
or the peak luminosity of SAX J1808.4-3658, but may simply indicate a low mass accretion rate. 

A definitive association of SAX J0635-I-0533 with the EGRET source requires detection of a 
periodicity in gamma-rays at the pulsar spin period. Due to the long integration time required to 
obtain a detectable gamma-ray signal, only a priori knowledge of the binary parameters would 
permit a sensitive search for periodicity in gamma-rays. This can be obtained with additional 
X-ray observations of SAX J0635+0533. 



We wish to thank Enrico Massaro (University of Rome) and Ignacio Neguerela (BeppoSAX 
SDC) for scientific discussions, and Guido Vizzini (IFCAI) for the technical support in the data 
reduction procedure. 
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Fig. 1.— Finding chart for SAX J0635+0533 from the Digitized Sky Survey in a field of 4' x 4'. 
The X-ray error circle has a radius of 30" (see text). 

Fig. 2.— a) The MFCS light curve of SAX J0635+0533. b) Hardness ratio vs time between 4- 10 
and 1.8-4 keV energy bands. 

Fig. 3.— The Z^(u = 26) plot for MECS+LECS events. 

Fig. 4. — The chance occurrence probability to find values greater than Z!f(i' = 26) (continuous 
line) compared with the measured integral distribution (dots). 

Fig. 5. — Pulse profile obtained taking into account / and /. Two cycles are shown for clarity. The 
superimposed error bar correspond to 1 standard deviation. 



T — \ — I — I — \ — I — r 



— I — I — I — I — \ — I — I — I — I — \ — r 



o 

X 
CD 



o 

X 



CD 



O 

X 

cv 



\ \ I \ I I \ I I \ \ 

CO C\2 ^ 

odd 

s /s^uno3 



_U I I I I \ I I I I \ L 



CO C\} 



o 



0^ 



ID 



PL, 



d 



J I I I I I I I I I I I I I I I I U I I L 

009 OQt OOt 09G 



